Noise Correlation Scalings: Revisiting the Quantum Phase Transition in 
Incommensurate Lattices with Hard-Core Bosons 
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Finite size scalings of the momentum distribution and noise correlations are performed to study 
Mott insulator, Bose glass, and superfluid quantum phases in hard-core bosons (HCBs) subjected 
to quasi-periodic disorder. The exponents of the correlation functions at the Superfluid to Bose 
glass (SF-BG) transition are found to be approximately one half of the ones that characterizes the 
superfluid phase. The derivatives of the peak intensities of the correlation functions with respect 
to quasiperiodic disorder are shown to diverge at the SF-BG critical point. This behavior does not 
occur in the corresponding free fermion system, which also exhibits an Anderson-like transition at 
the same critical point, and thus provides a unique experimental tool to locate the phase transition 
in interacting bosonic systems. We also report on the absence of primary sublattice peaks in the 
momentum distribution of the superfluid phase for special fillings. 
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I. INTRODUCTION 

Since the celebrated work of Anderson [l| and others 
about half a century ago, the subject of disorder and 
quantum phase transitions induced by disorder contin- 
ues to attract interest Q ■ In spite of being intensively 
studied, the interplay between interactions and disorder 
remains an open frontier. In recent years, a new im- 
petus toward the understanding of disorder systems has 
emerged in the context of ultracold atomic systems be- 
cause of their extraordinary degree of tunability. In those 
systems, it is now possible to control the degree of dis- 
order either by imprinting speckle patterns 0, or by 
superimposing a secondary lattice on the main lattice 
to generate a quasiperiodic potential Using those 

techniques, Anderson localization has already been ob- 
served in one dimension At the same time, the 
possibility of exploring strongly interacting regimes with 
ultra-cold atoms has been demonstrated by observation 
of the superfluid to Mott insulator transition |9l-ll2j. and 
the realization of the Tonks-Girardeau gas, i.e., a gas of 
impenetrable (hard-core) bosons (HCBs), in one dimen- 
sion [l3l - [l5| . The combined investigation of the effect 
of disorder and correlations in low-dimensional systems 
is therefore a window that is now open for experimental 
exploration flfj| and has driven many recent investiga- 
tions [T3-H31. 

In this paper we revisit the problem of one-dimensional 
(ID) HCBs with pseudo-random disorder generated by 
imposing two-color lattices with periods that are mu- 
tually incommensurate. Even in the absence of inter- 
actions, one-dimensional quasiperiodic systems exhibit a 
localization-delocalization transition. In contrast to truly 
random ID systems where states are localized irrespec- 
tive of the strength of disorder, quasiperiodic systems ex- 
hibit a quantum phase transition from extended to local- 



ized at a finite strength of the quasiperiodic disorder. In 
the interacting soft-core regime, several studies have ad- 
dressed the phases that result from the interplay between 
pseudorandom disorder and strong correlations. Those 
studies have utilized different com put ational techniques 
such as exact diagonalization [2ll. |22|| . quantum Monte 
Carlo simulations [23[ , and the density- matrix renormal- 
ization group method 0, HU . A recent comprehensive 
discussion of these, and other, results can be found in 
Ref. fl6j . In the present work, we perform an exact nu- 
merical study of ultracold bosons in the hard-core regime, 
as well as spinless fermions, and investigate the effects of 
the interplay between disorder and interactions in the 
momentum distribution and noise correlations of those 
systems. These two observables can be experimentally 
probed through time-of-flight measurements (2(| H3] • 

Noise correlations have been examined before for sev- 
eral HCB systems, including homogeneous, period-2 su- 
perlattices, and truly random disordered systems [28l — 
|30(. However, previous studies of HCBs in incommensu- 
rate supcrlattices 17[ were limited to small system sizes, 
which prevented an accurate determination of the scaling 
of off-diagonal correlations at the superfluid to insulator 
transition. The present study is partially motivated by 
an alternative approach to computing four-point corre- 
lations (30[. This approach enables the investigation of 
larger systems and opens the possibility of a detailed scal- 
ing analysis and better characterization of the quantum 
phase transition. One of our key results here is the deter- 
mination of the exponents at the transition point. Fur- 
thermore, as one approaches the thermodynamic limit, 
roots of divergent behavior in the derivative of the cor- 
relation peak intensities become apparent at the HCB 
localization-delocalization phase transition. That singu- 
larity is absent in the free fermionic system. This is an- 
other important finding in this work. 
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The paper is structured as follows. Following the back- 
ground introduction, Sec. II Al describes the model and 
defines the correlation functions. The localization tran- 
sition is then reviewed in Sec. II B[ Section [XT] is devoted 
to the scaling study of momentum distribution and noise 
correlations, revealing the distinctive behavior at the crit- 
ical point. In Sec. IIIII the dependence of both quantities 
on the strength of disorder is explored and a new char- 
acterization of the phase transition is given. Section IIVI 
reports on the peculiar absence of correlation features for 
some fillings. In Sec. |V]we summarize our results. 



A. Model 

One-dimensional hard-core bosons in incommensurate 
superlattices, within the one-band approximation, are de- 
scribed by the following Hamiltonian: 



H 



HCB 



H.c. 



in 



where bj(bj) is the bosonic creation (annihilation) op- 
erator at the site j, satisfying the bosonic commuta- 
tion relations [Sj,Sj] = 5ij and the on-site constraints 

b\ 2 = bf = 0; n b j = bjbj is the bosonic particle number 
operator. Here t is the hopping energy between adjacent 
sites, which we set to one [t = 1) in our work, and Vj 
denotes the external potential on site j generated by the 
secondary lattice. Note that the primary lattice is the 
only one that is described within the tight binding ap- 
proximation. Vj, on the other hand, has the following 
form: 



Vj = 2Acos(2ttct7 + 5), 



(2) 



where the parameter A is proportional to the intensity 
of the lasers used to create the secondary lattice, a is 
the ratio between the wave vectors of the two lattices, 
and 5 is a generic phase factor. Here, we choose 6 to 
be zero without lost of generality. Finally, we note that 
periodic boundary conditions are considered throughout 
this work. 

Quasiperiodic lattices are characterized by an irra- 
tional a, which results from the incommensurate period- 
icity of the two superimposed lattices. In our studies, we 
choose a to be the reverse golden ratio, a = (V5— 1)/2 = 
0.61803.... This choice is motivated by the fact that the 
golden mean is considered to be the most irrational num- 
ber [l9l . l3"lj and has simple number theoretical properties. 
In our numerical calculations, a is approximated by the 
ratio of two consecutive Fibonacci numbers Fm-i/Fm, 
where the Fibonacci numbers are defined by the recur- 
sion relation F i+ \ = F t + Fj_i, with Fi = Fq = 1. 
Within the approximation followed in this manuscript, 
model ([TJ on a quasiperiodic supcrlatticc reduces to a 
model on a periodic superlattice whose properties are 



well known (a is rational). The latter nevertheless cap- 
tures the localization-delocalization transition exhibited 
by the quasi-periodic system [H, HH . 

We focus our analysis on two observables, the momen- 
tum distribution function n k and the noise correlations 
Akk' j defined as 



(3) 



Afefe- = {h k h k ,) - (n k )(n k >) - (n k ) (h-k> 



,nK — "k.k' 



where fc(fc') denote the momentum values, L the num- 
ber of lattice sites, a the lattice constant, K = 2ir/a the 
reciprocal lattice vector, and n an integer. The expres- 
sion for A kk i was derived in Ref. j28j . The linear term, 
proportional to (n k ), has its roots in the commutation 
relations of bosonic fields and the requirement of using 
normal ordered operators if one wants to map time-of- 
fiight observables to in- situ correlations evaluated before 
the expansion and restricted to the lowest band. Sec 



Ref. 341 for details. 



An important point to be kept in mind when evaluat- 
ing it is that, for hard-core bosons, only normal ordered 
expectation values can be computed using the mapping 
to a spin- 1/2 chain and then to nonintcracting fermions 
[30| . In contrast to earlier studies where correlations 
were calculated using Wick's theorem [28j | , here we follow 
an alternative procedure based on Slater determinants 
[30l l35l l36j. which, for completeness, is briefly reviewed 
in Appendix [Al 



B. Localization transition 

As mentioned before, two-color superlattices charac- 
terized by an irrational a, are intermediate between pe- 
riodic and fully random, and exhibit a metal-insulator 
transition in one dimension. The critical point X c = t 
describes the onset of localization. Below criticality 
(A < A c ), single-particle states are extended, namely, 
they are Bloch-like states. Above the critical point, 
single-particle states are exponentially localized with lo- 
calization length £ = ln(A) -1 [37j . 

In the extended phase, the energy spectrum consists of 
a set of bands (it becomes a Cantor set at the localization 
threshold). The gaps are located at k = ±na(mod l)w/a 
in the first Brillouin zone, indexed by an integer n. Note 
that as n increases the width of the gaps narrows. As 
a consequence of our selection of a — Fm-i/Fm, the 
single-particle spectrum consists of Fm bands and Fm — 1 
gaps. The filling factor p = N/L in a many-body system, 
where N is the number of particles and L the number of 
sites, is then a control parameter that tunes transitions 
between different phases. For fillings p = ncr(mod 1) 
and p = 1 — ncr(mod 1), the ground state of HCBs in 
the extended regime is a Mott insulator. This insulating 
state is incompressible and it is also referred to as an 
incommensurate insulator uj{ H9, 23] . The fact that the 
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FIG. 1: (Color online) (a) Momentum distribution n(k) and (b) noise correlations A^o, of HCBs, as a function of k for three 
different values of A. In all cases L = 233, iV = 166, and a = 144/233. The insets in (a) and (b) show the set of sublattice 
peaks in the momentum distribution and noise correlations, respectively. 



system is insulating at those fillings can be understood 
because the conduction bands in the fcrmionic system, to 
which the HCBs can be mapped, are filled [13, EH . For 
all other fillings, the ground state is supcrfluid. To be 
able to observe the fractional Mott insulating phases, the 
lattice size must be identical to (or multiple integers of) 
the period of the incommensurate potential, which means 
that we choose L = Fm while having a = Fm—i/ -Fjw • In 
our study, we consider lattices with L =34,55,89,144, and 
233. For A > A c , i.e., in the localized regime, the ground 
state remains an incompressible insulator for fillings p = 
ner(mod 1) and p = 1 — ncr(mod 1). For all other fillings, 
on the other hand, the ground state is a Bose glass. The 
Bose-glass phase is a compressible insulating state where 
the absence of transport is the result of localization [T|| . 

II. QUANTUM CORRELATIONS: SCALING 
ANALYSIS 

In this section, we perform a finite size scaling study 
of the momentum distribution as well as the noise corre- 
lations in different regimes. Typical momentum profiles 
and noise correlation patterns are shown in Fig. [TJ They 
were calculated on a system with 200 sites at half filling 
and for three different values of A representing the super- 
fluid phase, the critical point, and the localized phase. 

The momentum distribution reflects the quantum co- 
herence of the system. It exhibits a sharp diffraction 
pattern in the extended regime (supcrfluid phase) and a 
flat profile in the Bose glass and Mott insulating phases 
[ see Fig. QJa)]. The interference pattern is character- 
ized by a peak at k = and additional sublattice peaks 
at the corresponding reciprocal lattice vectors. In the 
quasiperiodic lattice under consideration, they are found 
at ka = ±n<7(mod l)27r and ka = ±(1 — ?ier)(mod l)27r. 
In general, only peaks associated with the main recipro- 
cal lattice vectors, n = 1, arc distinguishable. The k = 
peak is the sharpest peak and signals the supcrfluid char- 
acter of the extended phase. 

The noise correlation pattern contains information 



of the momentum-momentum correlations and is also 
known as Hanbury Brown-Twiss interferometry. It ex- 
hibits interference peaks in all phases, even in the insu- 
lating ones. In contrast to the momentum distribution, 
the noise correlation fringes have their root in the quan- 
tum statistics of the particles. For bosonic systems, the 
peaks are positive and are a manifestation of bunching. 
For fermionic systems, they are negative and signal anti- 
bunching [26|, [27j • 

The height of the noise correlation peaks and their 
underlying background significantly vary in the various 
quantum phases studied here. For example, a negative 
background is found in the superfiuid regime. Also, the 
sublattice peaks disappear in the Mott insulting phase 
which occur at unit filling. The insets in Figs. [Ha) an d 
[TJb) show the most noticeable sublattice peaks in our 
systems. 

In Fig. [2] we plot the central peak of the momentum 
distribution function rik=o ( a ) an d the noise correlation 
peak Aoo (b) as a function of the system size L for var- 
ious values of A. Note that the filling factor is fixed at 
p = 0.5 so that a superfiuid state is associated with the 
A < 1 cases. The distinction in the scaling behavior in 
the extended regime, at the transition point, and in the 
localized regime is apparent. One can see that both nk=o 
and Aoo exhibit power-law scaling in the superfiuid phase 
as well as at criticality, while in the Bose-glass phase they 
remain essentially unchanged with changing L. 

In Fig. [H[a), the four sets of values of nk=o vs L for 
A =0.5,0.6,0.8, and 0.9 exhibit exponents 0.503 ± 0.001, 
0.496 ± 0.002, 0.462 ± 0.007, and 0.42 ± 0.01, respec- 
tively, which are compatible with the well-known rik=o ~ 
\f~L scaling for HCBs in the supcrfluid regime, where 
the quasi-long-range order is present in the one-particle 
correlations. In Fig. &b), the exponents for Aoo in 
the extended regime are 0.903 ± 0.003, 0.878 ± 0.004, 
0.799 ± 0.007, and 0.72 ± 0.01, also for A =0.5, 0.6, 0.8, 
and 0.9, respectively. These results are compatible with 
the linear behavior Aoo ~ L reported in Ref. [3(|, where 
the supcrfluid phase of the homogeneous system and the 
period-2 lattice was studied. Overall, the results in the 
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FIG. 2: (Color online) Scaling of (a) rik=o and (b) the central 
peak Aoo in half-filled systems for various values of A. The 
inset in (a) depicts the scaling of nk=o for A = A c and up to 
larger system sizes. The black dotted lines are power-law fits 
to the data, from which the exponents (see text) indicate the 
distinction between the extended regime, the transition point, 
and the localized regime. 



extended regime are consistent with the fact that for 
A < 1 the system is a Luttinger liquid with K = 1 [l6| . 
On the other hand, for A > 1, Fig. [2] shows that both 
rik=Q and Aoo do not scale with system size, which is a 
demonstration of the presence of exponentially decaying 
correlations in the Bose-glass phase, and hence the lack 
of Luttinger liquid behavior. 

The behavior at the critical point is the most interest- 
ing one, as the system is not expected to be described 
by the Luttinger liquid theory. Still, we find a power-law 
scaling of both rik=o and Aoo with an exponent that is 
roughly one-half of the one in the Luttinger liquid regime. 
The inset in Fig. HJa) shows nk=o for A = 1 with larger 
system sizes than those depicted in the main panel and 
results for an additional filling p = 1/3. We find that 



22±0.01 



at p = 1/2 and nk=o 



-0.17±0.01 



at p = 1/3. In Fig. [UJb), A 00 at the transition point 
scales with an exponent of ~ 0.48 ± 0.03 at p = 1/2 and 
~ 0.45 ±0.02 at p = 1/3. Notice that there are apparent 
fluctuations in the data, at the critical point and above, 
with increasing system size. These are finite size effects 
whose main origin may be the fact that the filling factor, 
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FIG. 3: (Color online) Scaling of the sublattice peaks A% of 
half-filled systems for various A. These plots show that the 
sublattice peaks do not scale with system size. 



as well as a, fluctuate slightly between lattices with dif- 
ferent size. Unfortunately, the number of particles and 
the number of lattice sites cannot be accommodated to 
give the same exact filling factor with increasing system 
size. Still, our results are consistent with a halving of 
the exponents from their Luttinger liquid values. This is 
something that deserves further theoretical investigation. 

We carry out a similar study for the sublattice peaks. 
The L dependence of the primary sublattice peak (n = 1) 
is presented in Fig.[3]for the same values of A as in Fig. [21 
Notice that the sublattice peaks are immersed in a neg- 
ative background, as shown in the inset Fig. QJb), so we 
quantify their height by subtracting the background; that 
is, we define A; = A a ^ fi - ^A^ +aM + K^-6k,o), 
with 8k = 2-k/Lci. For systems at p = 0.5, the sublat- 
tice peaks do not scale with system size for any value of 
A. This finding is surprising for A < A c because, as fol- 
lows from the scaling of the Aoo peak, quasi-long-range 
order is present in the system in that regime. On the 
other hand, the absence of those peaks is expected in the 
glassy phase where correlations decay exponentially. Fur- 
ther analysis reveals that, within a range of filling factors 
from p = 1 — a to p = a, A* does not scale with system 
size in the superfluid phase. This will be discussed in 
detail in Sec. [EVl 



At the incommensurate fillings p = n<r(mod 1) and 
p = 1 — no (mod 1), a Mott insulating phase rather than 
a superfluid one is present for all A ^ 0. We then study 
the noise correlation peaks in systems with incommensu- 
rate filling p — a. In Fig.lH both Aoo and A* are plotted 
as a function of L for various values of A. The absence 
of scaling is apparent in both noise correlation peaks for 
all regimes. For those fillings, and the noise corre- 
lations do not distinguish the localization dclocalization 
transition that occurs at A r . 
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FIG. 4: (Color online) (a) Central noise correlation peak Aoo and (b) the sublattice peak A% as a function of L in systems at 
incommensurate filling p = a. 



III. QUANTUM PHASE TRANSITION 

In order to gain a better understanding of the disorder- 
induced phase transition and its signatures in our 
strongly interacting system, we study the behavior of the 
noise correlations in the vicinity of the critical point. Our 
study involves monitoring the peak intensities and their 
derivatives as the strength of the quasipcriodic order is 
tuned. 

Figures H[a) and [S£b) depict the results for the cen- 
tral peak and the primary sublattice peak, respectively, 
as a function of A in a half-filled system. With the in- 
crease of disorder, the height of the central peak (A 00 ) 
decreases and attains an asymptotic value Aoo = p(p+l) 
as A — > co [T3, H3] , while the height of the sublattice peak 
increases with disorder and saturates at a finite but small 
value as A —> oo. In the neighborhood of the the critical 
point, a kink is seen in both peaks as A is varied. In 
order to better understand the behavior in the critical 
region, we compute dAoo/dX. Figure [SJ^a) shows that 
the derivative of the height of the central peak exhibits 
a maximum at the transition point. A similar feature 
is seen in the derivative of the primary sublattice peak 
intensity [Fig. [SJb)]. 



It is instructive to compare the behavior of the noise 
correlation peaks at the localization transition in HCB 
systems with the one at the metal-insulator transition in 
the corresponding free fcrmionic system. Analogous to 
the HCB case, in an incommensurate lattice system, non- 
interacting spinless fcrmions exhibit a localization transi- 
tion at A c . Fingerprints of the phase transition in the cor- 
relation functions of the fermionic system have been sys- 
tematically studied in Ref. [l9j]. Our comparative study 
is aimed at singling out distinctive features in both sys- 
tems. Notice that the HCBs are interacting and have 
bosonic statistics, while their fermionic counterpart is 
nonintcracting and, of course, has fcrmionic statistics. 
Noise correlations are sensitive to both interactions and 
the statistics of the particles involved. 

Figure [5] shows the behavior of noise correlations in the 
fcrmionic system vs A. In this case, distinctive features 
of the noise correlations are the Bragg dip at ka = 2n 
[depicted in Fig. |U[a)] and the primary sublattice dip 
at ka = 2<77r [depicted in Fig. [Bfb)], as well as their 
derivatives with respect to A. As in the HCB case, the 
metal-insulator transition in the fermionic system is also 
signaled by peaks in the first-order derivatives of noise 
correlations at the critical point. 




FIG. 5: (Color online) (a) Central peak Aoo and (b) sublattice peak A* as a function of A and their corresponding derivatives 
in a system with N — 72, L = 144. The onset of the superfluid to Bose-glass transition is marked with a peak in the first-order 
derivative of both noise correlations. Note that in both panels the left scale applies to the intensity of the peak while the right 
one applies to its derivative. 
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FIG. 6: (Color online) (a) Primary dip \A2jr_ | and its derivative, and (b) the sublattice dip |Ao-| and its derivative, are plotted 
as a function of A for noninteracting fermions. N = 72, L = 144. In the fermionic system the critical point is also signaled by 
sharp peaks in the first-order derivatives. Note that again in both panels the left scale applies to the intensity of the dip while 
the right one applies to its derivative. 



A study of the scaling of the peaks and their derivatives 
with increasing system size makes evident that there is 
an important difference between their behavior for HCBs 
and fermions. The results of such a scaling, for the 
derivatives at the critical point, arc presented in Fig. [7] 
For HCB systems, we plot dA 00 (X c ) / dX and dA* a (X c ) / dX 
as a function of L. One can clearly see that dAoo(X c )/dX 
diverges with increasing system size. This is expected 
as, on the left of the transition point, Aoo increases with 
system size while, on the right of the transition, Aoo 
saturates to a A-dependcnt (size-independent) value as 
the system size is increased. As also shown in the fig- 
ure, the behavior of dn,k=o{X c ) / dX is qualitatively sim- 
ilar; dA*(X c )/dX, on the other hand, does not change 
with system size, as follows from the lack of scaling of this 
sublattice peak discussed with Fig. [3] For the fermions, 
shown in the inset in Fig. [7J one can see that none of the 
derivatives scales with system size. 

Hence, in the thermodynamic limit, the derivatives 
of both the HCB and fermion noise correlations peak 
heights exhibit maxima at the critical point. However, 
it is only in the HCB system that these maxima diverge. 
Hence, they provide very sharp signatures to locate the 
phase transition in experiments with ultracold gases in 
the bosonic case. 



IV. ABSENCE OF PRIMARY SUBLATTICE 
PEAKS 

In this section, we analyze the absence of scaling of the 
primary sublattice peaks, in the supcrfhhd phase, for a 
specific window of filling factors. 

As already discussed in the previous sections, there 
is an unexpected lack of scaling with system size of the 
sublattice peak A* when p = 0.5 in superfhhd phase. 
In Fig. [8] we plot the sublattice peaks of the momentum 
distribution and the noise correlations as a function of 
filling factor p for A = 0.5, i.e., within the supcrfluid 
phase. We show results for two system sizes to illustrate 



the scaling behavior with system size. Similar to A*, 
we define the primary sublattice peak in the momentum 
distribution as n* = — ^(n (J 2 2L+Sk + n (r zn__g k ), with 
5k = 2n/La. It is apparent in Fig. [SJa) that the sublat- 
tice peaks n* arc absent in the window of filling factors 
between the fractional Mott states at incommensurate 
fillings p = 1 — a and p = a, respectively. Correspond- 
ingly, for A* [Fig. Efb)], a system size invariance is ob- 
served within the same filling window. The absence of 
scaling in the noise correlations in that window of fillings 
then follows the absence of sublattice peaks in momen- 
tum distribution. 

As mentioned in Sec. II B[ for A < A c , the one-particle 
energy spectrum is split up by gaps at incommensurate 
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FIG. 7: (Color online) (Main panel) Derivatives of Aoo, A*, 
and nfc = o as a function of the system size at A c in HCB sys- 
tems. A divergence in both dAoo(A c )/dA and rik=o(Xc)/dX 
is seen as the system size increases, while dA* a (\ c ) / dX satu- 
rates to a constant value. The inset shows IdA^ /dX\ and 

\dA a /dX\ as a function of L, and also at A c , in fermionic sys- 
tems. No size dependence is seen in this case. The mean 
occupation per lattice site is p — 0.5 for both cases. 
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FIG. 8: (Color online) The primary sublattice peaks as a 
function of p in (a) momentum distributions and (b) noise 
correlations, for two different system sizes, A = 0.5. Within 
the range from pi = 1 — a to p2 = a (marked by blue dotted 
lines), the vanishing of peaks and the absence of scaling are 
noticed in n* a and A£, respectively. 



values of k and the two largest gaps at k = ±an/a are 
associated with the fillings p = 1 — a and p — a, re- 
spectively. As a result, the fillings within the window 
p = 1 — a and p = a correspond to supcrfluid states with 
a partially-filled central band which is sandwiched by the 
two a gaps. (There is actually a set of gaps in the cen- 
tral band, simply ignored as they are small compared to 
the a gaps.) Interference effects in this window of fillings 
then seem to prevent those sublattice peaks from emerg- 
ing in n(k) and from scaling with system size in the noise 
correlations. 

We should mention that a similar behavior can be 
seen in commensurate superlattices, but only in the case 
where the period is a multiple of 4. Interestingly, for 



those periods, it was discussed in Ref. 33[ that the two 
central bands in the spectrum cross at k = 0. Results 
for the momentum and noise correlation peaks for 
commensurate superlattices with period-4 and period-8 
are presented in Fig. [9] There, we plot Tip cl function 
of p for a = 1/4 and a = 3/8, respectively. The primary 
sublattice peaks in momentum distributions are absent 
at fillings between 0.25 and 0.75 for a = 1/4 [Fig. Ha)] 
and between 0.375 and 0.625 for a = 3/8 [Fig. |^b)]. 
Again, we attribute the missing peaks to a destructive 
interference of the single particle wave functions, which 
may be related to the central band overlaps. 

V. CONCLUSIONS 



Critical scaling exponents characterizing the behavior 
of correlation functions are the hallmark of critical phases 
and phase transitions. In this paper, we have character- 
ized the disorder-induced quantum phase transition be- 
tween a supcrfluid and a Bose-glass phase in two-color 
incommensurate superlattices by computing scaling ex- 
ponents of the momentum distribution and noise corre- 
lation peaks at k = 0. At the quantum phase transition, 
those obscrvables, which are accessible in cold gases ex- 
periments, are found to scale with system size as ~ TV 25 
and ~ A 5 , respectively. Intriguingly, the values of the 
exponents are one-half of their corresponding values in 
the superfluid phase, where the system is known to be 
described by the Luttinger liquid theory with Luttingcr 
parameter K = 1 [l6jj. We also show that, as expected, 
the scaling of those obscrvables with system size vanishes 
in the localized phase. 

We have studied the behavior of momentum and noise 
correlation peaks, as well as their first derivatives, as the 
superlattice strength is varied. We have shown that both 
in HCB and noninteracting spinless fermion systems the 
supcrfluid-Bose-glass transition is signaled by peaks in 
the first derivatives of the noise correlation. However, 
only in the HCB case do those derivative peaks diverge 
with system size, providing a sharp experimental signa- 
ture of the phase transition. 



0.20 



0.15 



0.05 



0.00 



A 


— •— L=80 
— ■— L=160 


1 1 

A 


if 'J 
f? V 


a=1/4 


■/v ■ 

If \ i : 


■ 

• 






■ (a) 




■ 







0.0 



0.2 



0.4 0.6 
P 



0.8 



1.0 




FIG. 9: (Color online) as a function of filling factor in commensurate systems with period 4-multiples: (a) a = 1/4 in a 
period-4 superlattice and (b) a = 3/8 in a period-8 superlattice. Two different sizes are plotted in both cases, for A = 0.5, 
which makes evident that the sublattice peaks are absent for fillings between p — a and p = 1 — cr. 
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Appendix A: Exact computational approach 

In order to solve the HCB problem, we map the Hamil- 
tonian described in Eq. ([1} to the exactly solvable non- 
interacting fcrmion Hamiltonian 



fi fi+1 



H.c. 



(Al) 



by means of a two-step process. The first step is to 
map HCBs to spin- 1/2 systems through the Holstein- 
Primakoff transformation (HPT) [H| 



1 - bib. 



1 



PA h 



(A2) 



where af 1 are the spin raising and lowering operators for 
spin-1/2 systems. The HPT leads to a replacement of 
&j(Sj) by o^~(ct7") if and only if the HCB creation and an- 
nihilation operators are normal ordered before the map- 
ping. This is a subtle point that needs to be considered 
for a correct calculation of the noise correlations if one 
wants the HCB Hamiltonian to be the limit U —> oo of 
the Bose-Hubbard model. 

Subsequently, the spin-1/2 Hamiltonian is mapped 
to the fermionic Hamiltonian by means of the Jordan- 
Wigner transformation (JWT) [39| 

*t = fff[e- iv ^, K = f[e™ fV *k (A3) 

13=1 (3=1 

in which fj (/j) is the creation(annihilation) operators 
for spinless fermions. In the nonintcracting fermionic sys- 
tem the ground-state wave function can be expressed as 
a Slater determinant 



N L 

K — l fJ=l 



(A4) 



with the matrix (P)l,n given by the lowest N single- 
particle eigenfunctions of the Hamiltonian in Eq. (|A1|) . 
From the fermionic Slater determinants one can then cal- 
culate the two spin-1/2 Green's functions, Gij and Gijki- 



G, 



i-l 3-1 

ft f_ 2 ;+ "n- _„•„ ft , 



13=1 7=1 



det 



P"F 



(A5) 



G 



ijkl 



a=l 13 = 1 

fc-1 (-1 

7=1 5=1 

det P w ] , (A6) 



which are required to compute rife and Akk' in Eq. Q 
[30j . In the expressions above, 



pa 



-Pg K for g < a, k = 1, . . . , N 
P gK for g > a, k = 1, . . . , N 
5 ae for k = N + 1 



(A7) 



and 



pap 



— P£, for < Ot, K = 1, 
1, 

7V + 2 



P£. for p> a, k 



..,N + 1 

,.,N + 1 (A8) 



S a g for K 



with a(/3) = i,j,k,l. The most time-consuming part of 
our calculations is determining all the nonzero elements 
of Gijki-, each of which involves the multiplication of a 
(N + 2) x L matrix and a L x (N + 2) matrix, which 
scales as (N + 2) 2 L, and then computing the determinant 
of the resulting (N + 2) x (N + 2) matrix, which scales 
as (iV + 2) 3 . 
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